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The VLF and ELF modes excited by a horizontal magnetic dipole (vertical loop) in the 
spherical shell between a finitely conducting earth and an isotropic sharply bounded iono- 
sphere are shown to have a nearly transverse magnetic character. The modes are similar 
to those of a vertical electric dipole. With the exception of the zero order mode, the propa- 
gating modes excited by the magnetic dipole are of slightly higher amplitudes, provided 
that the far fields of the horizontal magnetic and vertical electric dipoles are equal over flat 
earth in the absence of ionosphere. 

The transient fields generated by a current step in the magnetic dipole are in the first 
approximation similar to the fields generated by a current impulse in a vertical electric 
dipole. Response of the zero order mode of the magnetic dipole has been calculated. 



1. Introduction 

The mode theory of VLF traiisniissions has been 
developed by Watson [1919], Budden [1953], Sehii- 
mann [1954], and Wait [1957, 1960a, 1960b] for 
excitation by vertical and horizontal electric dipoles 
and also by vertical magnetic dipoles (horizontal 
loops). Additional references on this subject have 
been listed by Wait [1960a, 1960bl. 

Tlie reciprocity theorem lias been found to be use- 
ful for relating the fields of a horizontal electric dipole 
to fields of vertical electric and magnetic dipoles 
[Wait, 1960a]. The reciprocity relations nniy be 
also applied to relate the fields due to horizontal 
magnetic dipole excitation to known fields due to 
vertical electric and magnetic dipoles. 

2. Reciprocity Relations 

The reciprocity theorem states that the voltage V2 
induced in antenna 2 by current /i of antenna 1 is 
the same as the voltage Vi induced in antenna 1 by 
an identical current I2 flowing in antenna 2. 

A vertical electric (VE) dipole of length ds at ZQ=Zt 
and a horizontal magnetic (HM) dipole parallel to 
the a;-axis (or a vertical loop of area da in the y, z 
plane) at z^^'Zr are shown in figure 1, where sub- 
scripts T and i refer to receiver and transmitter 
coordinates respectively. The only nonzero H 
component generated by the VE dipole is H^'. The 
magnitude of the voltage induced in the HM. dipole 
yhm jg maximum if Il^> is parallel to the a:-axis, 
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which occurs for 0=^ and — • For e^'^^ variation of 
the fields. 

ynm^ -i^^^Hl'.(zo=Zt, Zo= Zn) sin da, (1) 

where the superscripts hm and ve refer to horizontal 
magnetic and vertical electric dipoles respectively. 
With the same current applied to the HM dipole 
only tlie vertical electric field E^ will contribute to 

V and 



y.e^JE^m^^_^^^^^^^^^^^^^ 



(2) 




♦• y 



FiGUBE 1. Coordinate systems for defining reciprocity relations. 
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Equating (1) {uul (2) results in 

= — ico/xo sin (l)Hl'f(zQ= Zi,ZQ= Zr)da/ds. (3) 

A vertical magnetic (VM) dipole (or a horizontal 
loop of area da^"^ parallel to the x^, y' plane) at 
z'Q = Zt and a horizontal magnetic {HM) dipole par- 
allel to the a:-axis (or a vertical loop of area da^"^ in 
the y, z plane) at z^^Zr may be considered in the 
geometry shown in figure 1. The only nonzero H 
component generated by the VM dipole is Hp' . The 
magnitude of the voltage induced in the HM dipole 
yhm jg maximum, if Hp' is parallel to the x-axis, 
which occurs for 0=0 or tt. Hence 



V'''^ = i0JfX,Hl'?{z'f,= 2t, Zo=Zr) QO^UoI" 



(4) 



where the superscript vm refers to the vertical mag- 
netic dipole. With the same current applied to the 
HM dipole only the vertical magnetic field H^ will 
contribute to V^^ and 



yv: 



-ioiixJH'l'^{2 



t,Zo=Zr)da''' 



(5) 



Equating (4) and (5) gives for da^'^ = da^ 



m^{z,= 2,,z',= Zr) = -^^OS<i>H;^{z,= Z,,Z,= Zr). (6) 

For spherical coordinates the z component is replaced 
by r component and the p component by the com- 
ponent of the fields. 

The reciprocity relations apply to dipoles in the 
presence of any linear media. The fields of dipoles 
above plane earth derived by solving the wave equa- 
tion should also satisfy (3) and (6). Thus (27) and 
(116) of Norton [1937], which have been derived for 
unit dipoles, check the validity of (3). Some algebra 
is involved in verifving (6) from (35.1) of Sommerfeld 
[1949] and from (2) to (6), (17) to (19), and (30) of 
Norton [1937]. 

3. Fields of the Horizontal Magnetic Dipole 

The dipole fields will be examined with the aid of 
reciprocity relations in an idealized spherical shell 
that is bounded by a homogeneous earth of radius a, 
conductivity a^ and dielectric constant eg and by a 
homogeneous ionosphere of radius (a+A), conductiv- 
ity (Ji and dielectric constant e^. Ihe above model 
neglects the effects of the earth's magnetic field. 

The vertical electric field El^'"" is computed by (3) 
applying (5.1), (6.16) of Wait [1960a] and using the 
relation 



de 



P,(— cos e) 



=-v[ms dP,(-cos ^)+P,_i(-cos ^)]/s'n e. (7) 



Tliis results in 
,f^^_ I(da)io)iJio sin </> 



E^' 



2hrr 



S /„,(.,)/„(..) -^^^-i^ 



[cos ^P,(-cos 6')+P,_i(-cos d)]. (cS) 



The vertical magnetic field Hi"' is computed bv 
(6) applying (9.1), (9.2) of Wait [1960a] and using 
the relation (7) as 



Hi 



I(da)oJfX0 cos </) ^ dmjm(^r)^ 

m=i sin MTTsin d 



2khrtrrif] 



\^m{^^)+r^J^J^{z^)\^ 



[cos dP,(~-cos e)+P,^i(-cos d)]. (9) 

Using the second order representation of the radial 
functions, applvhig (6.9), (6.17), (9.3), and (9.4) of 
Wait [1960a] yields 



2Uz)=e^'<^.^[Rl(C,)]-' 



/2_|.^-aC,.[^|(;Q]l/ 






(10) 



(11) 



where q is an integer equal to either norm and where 

(12) 



^„^^j_w/LJi 






^ m~\ '^p^p 



<7,=Vl-(sw^ 

p+0.5 



So 



ka 



(13) 

(14) 

(15) 

(16) 



The plus sign of (11) should be used with (Z=7i, 
while the minus sign is appropriate for q=m. In 
(16) p is a complex number, which is equal to v of 
(8) or equal to /x of (9), k=o)/c, c^velocity of light. 
The refractive index n of the boundary medium 
characterized by subscript p is 



:(o-^ + icoe^)/€o 



(17) 



where ap and e^ are the conductivity and permit- 
tivity of the boundary medium and where €o is the 
permittivity of free space (eo=(367r)"^ 10"^ farad/ 
meter). The boundary medium can be either ground 
(p=g) or ionosphere {p=i). The refractive index 
of ground is 



^(7^/(lcoeo) 



(18) 
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where a^ is llu* o-round ('oiiductivily. The refractive 
index of t he ionosphere is 

n\=l+a,l{iwe,). (19) 

TIk^ ionospheric conductivity a-, is defined by 

(20) 



enco^ 



where v is the electron collision frequency in the 
ionosphere and co^, is the plasma frequency (aj| = 3180 
A^, AT'^ number of electrons per meter^). Cn and 
Cm are roots of the modal equations 



m{C\)Rl{C\) = ex^ {+2ikh(\). 



(21) 



Sg that is related to Cq by (14) has a magnitude of 
jipproximately unity and^ a snudl imaginary part 
for propagating modes of low attenuation. With 
Sn iind S,„ (k^ter mined from (14), v of ((S) and /x of 
(9) follow from (16). This completes the formal 
specification of tlie fiehls E^T and //'';". 

The expressions iov E''^ and /Z^';" may be siin|)lified 
by introducing further approxinuitions. The con- 
sideration will be restricted to cases where 



(22) 

(2:]) 



This gives 



Vr-rt 


«fl., 


kC\Zr 


«1, 


kC,Z 


«1, 


ng> 1 . 


^2(a)«l- 


2 


?(a)«-l 




f /- -^ ~ 1 _i_ 


1 



fm{z) ~("'m/{in/\) -CJiin,,), 
(lf„(z)/dz=ikn/\f„Xz)~ikn,f„{z). 



(24) 
(25) 

(26) 
(27) 

(2S) 

(29) 
(30) 



Approximating the Legendre functions by [Watson, 
1919; Bremmer, 1949] 

Pp( — cos6/) «(27rp sind)-^-^ exp[i(p+~\7r-d) — iT/4l 

(31) 
results in 

£'r«=-shi cl>Et J^^ 
V sm (c( 






i/a) h/\ 



t:8„S',i'e ' ^ ^\ (32) 

71 = 



H^m^-COS 



.Eo"^ I d/a -\'rt 
rjUg V sin (d/a) h/ 



7/. . 27rrf . TT 



m = l 

where the distance along the curved earth is 
d—ad, 



and where 



^0 - ^2,; ^ 



(33) 

(34) 
(35) 



is the vertical electric field of the source at a distance 
d in the direction of maxinnnn intensity on a per- 
rectly conducting plane earth. 

4. Comparison of Harmonically Excited 
Dipole Fields 

The fields of the horizontal magnetic dipole Er"^ 
and //?'" will be compared first with the cori'espond- 
ing field components of the horizontal electric dipole. 
It follows from (9.37) and (9.42) of Wait [1960a] 
that for tlie n"^ and tlie m^^ mode 



fir 

TPhm 



^—Tig tan </) 



me 

J^ r 



=n, cot </> -~ 



(36) 



(37) 



where the supersci'ipts he and ve (h^sigmite fic^ld com- 
poiuMits of horizontal electric and of vertical electric 
dipoles respectively, /s'o"' is defined in (35), E'i^ is the 
corres])on(ling expression of the vei'tical electric (h])ole 
(h^fined by (6.25) or (9.3S) of Wait |1960a]. 1'he 
appi'oximation ^^^=^1 was used in the expressions 
for E^'/ and i??^ in order to make them consistent with 
the derivation leading to (32) and (33). Q~l con- 
stitutes a better approximation than/„(3)?^l and 
the assumption jn{^)^\ [Wait, 1960al should be 
folio w^ed by Cg^l. 

Comparison of (32) and (33) shows that r?//^'^ is 
several orders of magnitude snudler than A^?'"^ because 



and because 



n-^<l 



(38) 
(39) 



for propagating modes of low attenuation. The 
H^^ fields of a given mode are proportional to 
n~^-^-\l(jolaf,. The fields are decreased with decreas- 
ing frequency co or increasing ground conductivity 
ag. wSimilar comments apply to the transverse electric 
(TE) field components that may be derived from 
Hr. The TE fields will not be considered in more 
detail. 

The component Er"^ that is associated with trans- 
verse magnetic (TM) fields may be compared with 
the corresponding component Er^ of a vertical elec- 
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trie dipole using an expression of the Er^ fields de- 
rived under similar assumptions [Wait, 1957]. Thus 



Er 



sin (/)£'?^ 



ET 



~'SJ''ds 



sin <t>EV 



(40) 



Provided that the fields of the two antennas are the 
same in the given direction over a perfectly conduct- 
ing plane earth (sin </)£'J'^=£'g^), (40) may be simpli- 
fied to 



Ef" 



— ~~Sn- 



(41) 



Substituting the reciprocity relation (3) in (40) fol- 
lows that {El')l(y]H;') also satisfy (41). The latter 
ratio, obtained from (13), (15), and (16) of Wait 
[1960b], when used in conjunction with (3) provides 
an alternate way of obtaining (40) and subsequently 
(41). 

A more exact treatment of the horizontal magnetic 
dipole fields cannot be expected to lessen the im- 
portance of TM relative to TE field components. 
A better approximation to the TM fields of the hori- 
zontal magnetic dipole may be worked out from the 
exact ratio E^^/Er'"' derived in the appendix. As long 
as the source and observation point are at equal 
radii r the ratio involves only the azimuthal 6 func- 
tions explicitly. However, the accuracy of the roots 
Sn and of E^^ depends on the approximations to the 
azimuthal and radial functions. 

Numerical values of Sn are available from mode 
calculations of the vertical electric dipole [Howe 
and Wait, 1957]. For 0-^=00 the values are as 
follows: 



TTJ^CKcolh) 


3.5 
0.1 


3.5 
0.01 


7 
0.1 


1 

4 


0.9975-2 3X10-4 
. 98-2 3X10-3 
. 93-2 0. 01 
. 87-2 0. 02 


0. 997-i 10-8 
.97-28X10-8 
.91-2 9X10-8 
. 85-2 0. 01 


0.995-2 6X10-5 
. 992-2 4X10-1 
. 98-2 10-8 
.985-2 2X10-8 



The n=^0 mode is severely attenuated for fre- 
quencies above a few kilocycles [Wait, 1957]. In 
the lower frequency range the first order perturba- 
tion solution of the modal equation may be used 
[Schumann, 1954], which gives 



^n 



+AS:Li V^~^"+iV ^~^^J ■ ^^^^ 



The cutoff frequency of the nth mode is 



nwc 
''~h' 



(43) 



The perturbation method applies iovn9^0 if Cn^cjojco. 
Substituting this for Cn^ in the right-hand side of 



(42) and expanding the square root expressions 
,00 hoinLx^h ^i/ 2V< nyj 



^n 



1- 



IC 



2hc, (nl'^-' 



4) 



(44) 



The latter expression may be also obtained by solving 
(42) as a quadratic and by subsequently ignoring the 
second order perturbation. Applying (14), (42) 



simplifies for n=0 and 



\ni\:>l to 



^o=-\/l+T^ [-+-]• (45) 



For ng>ni and a-, = 1.2X10" 
has the following values 



10* 



10^ 



\ccllv = 1.35Xm So 



W 



1.04 6-^"2° 1.14 6-^70 1.46e-^'6« 



The above tabulations show that \Sn\ or the ratio 
lEr^jEf'^l are decreased with increasing mode index 
n and, with the exception of Sq, also with decreasing 
frequency. The phase angle between the propa- 
gating modes of the two dipoles is negligible. The 
largest difference between the fields of the two dipoles 
may be expected for the lowest propagating fre- 
quencies or near the mode cutoff ol n^O modes. 
The values of Sn near the mode cutoff may be deter- 
mined from an approximate solution of the modal 
equation. The first order perturbation solution is 
applicable only to the lower order modes. A differ- 
ent approximate solution is obtained by observing 
that Sn is small near the cutoff {Sn=l/{2 ka) at the 
mode cutoff for perfectly conducting ground and 
ionosphere according to Schumann [1954]). Hence 
one may look for solutions of the modal equation of 
the form 



Cn=^l-S'n=l-A, 



(46) 



where |AJ<1. Substituting (46) in (12), letting 



p=i and defining 



1=^^-^=% 



(47) 



results after neglecting terms with A^ and higher 
powers of A„ in 

^?=i(2i-i-2VlVi=^(l-^^5)=ii+A„Z2. 

(48) 

Considering 7i|»l, substituting (46) in (12) and 
letting p=g gives 

Rl^l-2n-'e'^''{l+A,)=Gi+AnG2. (49) 
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The solution of the modal equation (21) becomes 
2i{kh—'im) — InLiGi 



A„=- 



M Cri 



(50) 



It is simplified at the cutoff of the respective mode 
by applying (43) to 



An=- 



InLiGi 



2mx+|j+| 



(51) 



The calculated curves of \Sn\ and of its argument are 
depicted in figures 2 and 3 for infinite ground con- 
ductivity. \Sn\ is smallest for small values of 
I=Trcv/{hool) (or for large ionospheric conductivities) 
and for the higher modes (higher cutoff frequencies 
o^n^^nTrcjh). The argument of Sn exhibits a rela- 
tively small variation with n. The increase of [S'^l 
and of arg Sn with increasing / signifies the increase 
of mode attenuation with decreasing ionospheric 
conductivity. The increase of \Sn\ from the (Tg= oo 
to the (rg=0.07/hjcm curve (o-^=10~^ mho/m for h=70 
km) is less than 20 percent for the n=l mode in 
figure 4. This relative increase is even smaller for 
the higher n values. 

The amplitude ratio between electric and magnetic 
dipole fields in (41) will be larger for propagating 
modes than near the mode cutoff. The cutoff ratio 
is larger than 0.25 for the n=l mode and larger than 




Figure 2. Magnitude of the constant Sn near mode cutoff for 
perfectly conducting ground. 



0.15 for the first 8 modes with />0.01 (or for 
o■^> 1.2X10"^ with h=70 km). It may be noted 
that this comparison is made for dipoles that exhibit 
equal far fields above a perfectly conducting plane 
ground in the absence of ionosphere. The knowledge 
of the above amplitude ratios would enter in a com- 
parison of antenna radiation efficiencies. However, 
no such calculations have been attempted in this 
paper. 




Figure 3. Cutoff arguments of the! constant S^ for perfectly 
conducting ground. 




Figure 4. Magnitude of the constant S^ near cutoff for the 
n=l mode. 
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5. Transient Fields 

The transient fields of the horizontal magnetic 
dipole may be computed as the Fourier integral of 
(32) for a specified spectrum of the antenna current. 
For a step of current the far field approximation (32) 
gives 

E:^^{o,)=K^i^JZKSl-' e'"'^^'' (52) 



with 



i^— 



sin (f) Iq dar] 



(53) 



2hc^'27rac sin (d/a) 
and with 5o=0.5 and dn^o=^' 

The response of the n-=0 mode may be reduced 
to tabulated Fourier integrals after approximating 
S^o of (52) by a power series. With ?7^»?2i»l (45) 
reduces to 



So-^H 



h-\!io)iJi(Ti 



(54) 



Expanding VS'o and So in powers of (l/V'^<^) and 
ending the expansions with the quadratic terms 
results in 

ia)d ^liud 



Er^'^{o^)=Ke 



d/(2r/S2) 



V^*^" 



1 



2/3 8i82V"^w, 



;] 



where 



fi=2h^lX(Ti. 



(55) 
(56) 



Applying the transform pairs (806), (801), and 
(807) of Campbell and Foster [1948] gives 



7^,i/M-o.25/i r I I / 1\ 3"1 



where 



a;=T/a, 
T = {t-dlc)>0, 

e^ vd} 



«=(^ j =4t; 



A^ 4wW 



2C/32 8. 



d 



d ^J e; 



(58) 
(59) 

(60) 
(61) 



The power series expansion of (54) is applicable to 
frequencies where a)>coo=(^Vo'0~^. Hence (57) will 
be most accurate for T<;<^27r/coo or for 



a;<<Xo==27r 



[^va' 



=TrU, 



(62) 



Considering h^90 km, ^/=^10^ km, (7,==10~^ mho/m 
as an example, Xo=2.3 with 7i = 2.4. The representa- 
tion of the transient peak in figiu'e 5 will be accurate 
even for this pessimistically low a,-, value at d=10'^ 




Figure 5. Response to a current step of the n = mode of a 
horizontal magnetic dipole. 



km. Higher o-, and smaller d result in larger Xq 
(and u), which improves the accuracy of the transient 
tail representations. 

It follows from reciprocity relation (3) that for 
equal current waveforms P^{t) and I^^(t) in the 
HM and VE dipoles 



E^^{t) =—jjiQHl'{t) sin ct^da/ds. 



(63) 



A current in the VE dipole that is proportional to 
the time derivative of the current in the HM dipole 



results in 



E'{t)=TP'^{t), (64) 



£'^(0=--Mo sin ct^daHl'(t)l(Tds). (65) 



The response Er'^(t) to a current step Iou{t) should 
be proportional to the response H^it) to a current 
impulse IoT8{t). A^ssuming that (64) applies and 
assuming further that 



S,^l 



(66) 



(this is obviously incorrect at the lower frequencies), 
it follows from (40) that 



£"?"■(<)« -sin 4>Erit)daKcTds). 



(67) 



The impulse response of a vertical electric dipole 
has been calculated considering factors proportional 
to finite powers of So equal to unity [Schumann, 
1952 ; Wait, 1960bj. The calculated impulse response 
of the VE dipole [Schumann, 1952; Wait, 1960bl 
will be compared with the step response of the HM 
dipole (57) by means of (65) and also by means of 
the approximate relation (67). The leading term 
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of the impulse response (42) of Wait [1960b] (eqs 
(38), (39), (40), and (42) of Wait [1960b] should be 
nuiltipliod by 7r/2), when substituted in either (65) or 
(67) results 'in E^'\t) that is equal to (57) with the 
two terms proportional to u"^ and u~'^ of the square 
brackets set equal to zero. There is no agreement 
betw(HMi the higher response terms because of the 
S() approximation as discussed earlicM*. The impulse 
response of Schumann [1952], when substituted in 
(67) is the same as (57) with u= oo . 

Tlie first order perturbation solution of the modal 
equation becomes inaccurate for the higher fre- 
quencies involved in calculating the transient re- 
sponse of the 719^^ modes. Even the simplest (and 
inaccurate) expressions of Sn that may be obtained 
from (44) and (14) for |%|>>|/i,|>>l involve 
integrals similar to those encountered in the transient 
amdysis of lossy rectangular waveguides [Cerillo, 
194S]. The transient response of the tif^O modes 
has not been calculated in this paper. 

6. Appendix. Radial Electric Fields of 
Vertical Electric and Horizontal Magnetic 
Dipoles 

The fields of a vertical dipole nniy be derived from 
a single scalar function. Applving (29) of Schumann 
[1954] to (3) 



El' 



ds,§-^[.i 






"^ rF cos (j) da. sin 6 -^ 



m 



where for Uo the coordinates of the source and of the 
observation points are interchanged relative to Ui. 
As long as the source and the observation point are 
approxinnitely at the same radius {r^^rr), both u 
functions are the same. The ratio of the vertical 
electric fields of horizontal magnetic and vertical 
electric dipoles for equal currents on a perfectly con- 
ducting plane ground is 



El' 



Er 



ids'] 
kduo 



(69) 



Witli 



(68) becomes 



EV 



El 



■'1\ (-cos e) (70) 



(jz+l) sin e P^ (—cos e) 



E^'" El'" cos 4> rk cos 6 1\ (- cos 9) +P,_ , (- cos d) 

(71) 

The earlier approximations will reduce (71) to (41). 
However, more accurate v values [Wait, 1960a, ch. 
12] and a better approximation to /\ (—cos d) may 
be used in (71) to obtain a jnore exact E^^jE^"^ ratio. 
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tion coefficients, J. R. Johler, L. C. Waiters, and J. D. 
Harper, Jr., NBS Tech. Note 69 {PB 16 1570), Uuly 1, 1960) 
$2.00. 

The results of extensive computations performed during the 
course of a theoretical investigation of a sharply bounded 
model ionosphere for low- and very low-radiofrequency wave 
propagation are presented in the form of graphs and tables. 
This Tech. Note supplements the work described in, On the 
theory of reflection of low- and very-low-radio frequency 
waves from the ionosphere, J. Research NBS 64D, No. 3, 
269 {May-June 1960).];, 

The relation of /i,„«^F2 to M(3000)F2 and h„F2, J. W. Wright 
and R. Vl McDuffie, J. Radio Uesearch L(d)or(Uories 7, No. 32, 
409 {July 1960). 

Several easy methods of inferring directly from ionograms 
the height of the F region peak, using the quantities 
M(3000)/^2 or h^F2, are compared with values of hmaxF2 
obtained from N {h) profiles. It is found that such methods 
are free from bias at night at low and medium latitudes, 
but that systematic differences are significant in the daytime, 
and at high latitudes at all times. Examples are given 
illustrating the incorrect diurnal variation of hmax^^ estimated 
from these simple parameters, in comparison with simul- 



A test of a procedure for easy estimation of representative 
monthly electron density profiles for the ionosphere, J. W. 

Wright, /. Geophys. Research 65, No. 10, 3215-3217 {Oct. 
1960). 

A recently proposed method for obtaining a representative 
monthly electron-density profile for a given hour of the day 
applies standard N{h) procedures to a single representative 
virtual-height curve obtained in a special way from the 
individual daily virtual-height curves for that hour. The 
method is tested by comparing examples of the resulting 
profile with the average of the profiles for the individual 
days. The differences are smaU compared with the disper- 
sion of the individual profiles about their mean. It is con- 
cluded that the method is valuable for obtaining a world- 
wide N{h) morphology with a minimum expenditure of effort. 

Spiral patterns in geophysics, V. Agy, /. Atmospheric and 
Terrest. Phys. 19, No. 2, 136-140 {Oct. 1960), 
Recent analyses of magnetic and ionospheric data, notably 
by A. P. Nikolski, have resulted in spiral ''precipitation" 
patterns which lead the authors to claim support for Stormer's 
theory of the aurora. Although it may be possible to argue 
against these claims by attacking the methods of analysis 
and/or the Stormer theory itself, an entirely different approach 
is used in the paper presented here: an examination of per- 
tinent points of Stormer's theory shows that the analytical 
spirals mentioned cannot be Stormer spirals and this con- 
clusion holds regardless of the soundness of the analyses and, 
indeed, of the validity of Stormer's theory. 

Supplementary world maps of F2 critical frequencies and 
maximum usable frequency factors, Donald H. Zacharisen, 
NB8 Tech. Note 2-2 {PB151361-2) {Oct. 1960) $3.60. 



This report supplements NBS Tech. Note 2 (April 1959), 
and completes the basic data required for r2-layer maximum 
usable frequency predictions. Prediction charts are given 
for the months of February, April, May, August, October^ 
and November. Auxiliary charts are included to aid in 
predicting F2-layer MUFs. 

The four parameters used for predicting MUFs are foF2 and 
the 4000 km MUF factor for a twelve-month running average 
Zurich sunspot number of 50, and the rates of change of foF2 
and 4000 km MUF factor with sunspot number. The first 
three parameters are presented in map form for each even 
hour of Greenwich Mean Time. The fourth parameter is 
presented on a chart of geomagnetic latitude and local time. 

The height of maximum luminosity in an auroral arc, F. E. 

Roach, J. G. Moore, E. C. Bruner, Jr., H. Cronin, and S. M. 
Silverman, J. Geophys. Research 65, No. 11, 3575-3580 
{Nov. 1960). 

The height of maximum luminosity of an auroral arc is 
estimated from simultaneous obscu'vations at three stations 
in western United States during a night of general auroral 
activity (November 27-28, 1959). Photometricahy this arc 
is characterized by a selective enhancement of the [01] 6300 A 
line. From twenty-four individual measurements the height 
is found to be 412 km with a standard deviation of ±23 km 
for one observation and ±5 km for the mean. The geo- 
graphical position of the arc, its orientation, and its movement 
during the night are discussed. 

FM and SSB radiotelephone tests on a VHF ionospheric 
scatter link during multipath conditions, J. W. Koch, W. B. 
Harding, and R. ,1. Jnnsen, PRE Trans. Commun. Systems 
CS-8, No. 3, 183-186 {Sept. 1960). 

Experiments have been carried out on an ionospheric-scatter 
link to observe the effects of long-delayed multipath signals, 
caused by F2 propagated back scatter, on the intelligibility 
of voice comnuniication. Frequency modulation and single- 
sideband modulation equipments were used for the tests. 
During periods when the back-scatter signal levels approached 
the level of the normal ionospheric-scatter signals, the 
frequency-modulation voice transmissions were unintelligible; 
however, under the same conditions, single-sideband voice 
communication intelhgibility remained at almost 100 per 
cent although there was some loss in quality. 

Radio refractometry, Jack W. Herbstreit, NBS Tech, Note 66 
{PB 161567) {July 1960) 50 cents. 

The optical refractive index is known to be determined 
principally by the temperature and pressure of the atmos- 
phere, whereas the radio refractive index is, in addition, 
affected by the water content of the atmosphere, the relation- 
ship between these quantities being expressed in the following 
way: 

N= (n-1) 10^=(^) (P+[48106/T]) (1) 



where total air pressure P, and water vapor e, are in millibars, 
and the temperature T is in degrees Kelvin [Smith, 1953]. 
The quantities P, e, and T have long been routinely measured 
at the surface of the earth by standard weather bureau sta- 
tions. For some time they have been measured from the 
surface up to great heights using balloon-borne radiosonde 
equipment at a large number of places over the earth's 
surface. More recently, equipment has been developed to 
measure rapidly and directly the radio refractive index of the 
atmosphere using radio techniques. The measurement of the 
radio refractive index properties of the atmosphere and the 
application to radio propagation problems is the subject of this- 
paper. 
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Journal of Research, Vol. 65 A, No. 2, March-April 1961. 
70 cents. 

Mass spectra of some deuteroethanes. Edith I. Quinn and 

Fred L. Mohler. 
Heats of hydrolysis and formation of potassium borohydride. 

Waher H. Johnson, Richard H. Schumm, Isa H. Wilson, 

and Edward J. Prosen. 
Heat of combustion of borazine B3N3H6. Marthada V. 

Kilday, Walter H. Johnson, and Edward J. Prosen. 
Thermodynamic properties of thorium dioxide from 298 to 

1,200 °K. Andrew C. Victor and Thomas B. Douglas. 
Calculated energy dissipation distribution in air by fast 

electrons from a gun source. John E. Crew. 
Vitrons as flow units in alkali silicate binary glasses. Leroy 

W. Tihon. 
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type AB4O9. Carl R. Robbins and Ernest M. Levin. 



Journal of Research, Vol. 65B, No. 1, January-March 1961. 
75 cents. 

On transient solutions of the '^baffled piston" problem. F. 

Oberhettinger. 
Special types of partitioned matrices. Emilie V. Haynsworth. 
Bound for the p-condition number of matrices with positive 

roots. Philip J. Davis, Emilie V. Haynsworth, and Marvin 

Marcus. 
Some computational problems involving integral matrices. 

Olga Taussky. 
Computational problems concerning the Hilbert matrix. 

John Todd. 
Index to the distributions of mathematical statistics. Frank 

A, Haight. 
Selected bibliography of statistical literature, 1930 to 1957: 

IV. Markov chains and stochastic processes. Lola S. 

Deming and D. Gupta. 
Building code requirements for reinforced masonry, NBS 

Handb. H74 (1960) 15 cents. 
Amplitude, probability distributions for atmospheric radio 

noise, W. Q. Crichlow, A. D. Spaulding, C. J. Roubique, 

and R. T. Disney, NBS Mono. 23 (1960) 20 cents. 
Standard frequencies and time signals from NBS stations 

WWV and WWVH, NBS Misc. Pub. 236 (1960) 10 cents. 
Research highlights of the National Bureau of Standards, 

Annual Report, fiscal year 1960, NBS Misc. Pub. 237 

(1960) 65 cents. 
Quarterly radio noise data — December, January, February 

1959-60, W. Q. Crichlow, R. D. Disney, and M. A. Jenkins, 

NBS TN18-5 (PB151377-5) (1960) $1.75. 
Quarterly radio noise data — March, April, May 1960, W. Q. 

Crichlow, R. D. Disney, and M. A. Jenkins, NBS TN18-6 

(PB151377-6) (1960) $1.75. 
On the nature of the crystal field approximation, H. Gold- 
berg and C. Herzfeld, NBS TN67 (PB161568) (1960) 

$2.50. 
Vapor pressures of organic compounds in the range below one 

milhmeter of mercury, E. E. Hughes and S. G. Lias, NBS 

TN70 (PB161571) (1960) 75 cents. 
Some experiments on the deposition of gasses at 4.2 °K, T. 

Baurer, NBS TN73 (PB161574) (1960) $1.00. 
Scattering of cobalt-60 gamma radiation in air ducts, C. 

Eisenhauer, NBS TN74 (PB161575) (1960) 75 cents. 
Soviet research in field electron and ion emission, 1955-59; 

an annotated bibliography, T. W. Marton and R. Klein, 

NBS TN75 (PB161576) (1960) $1.25. 
Oblique incidence receiving antenna array for a relative 

ionospheric opacity meter, A. C. Wilson, NBS TN78 

(PB161579) (1960) 50 cents. 



The dynamic compressibility of a rubber-sulfur vulcanizate 
and its relation to free volume, J. E. McKinney, H. V. 
Belcher, and R. S. Maivin, Trans. Soc. Rheology IV, 
347 (1960). 



Procedure for the determination of the noble metal content 
of dental gold alloys, H. J. Caul, W. S. Clabaugh and M. E. 
Susa, J. Am. Dental Assoc. 61, No. 4, 339 (Oct. 1960). 
Theory of anisotropic fluids, J. L. Ericksen, Trans. Soc. 

Rheology IV, 29 (1960). 
Rapid frequency analysis of fading radio signals, J. M. Watts 
and K. Davies, J.^ Geophvs. Research 65, No. 8, 2295 

(Aug. 1960). 
On Stokes flow about a torus, W. H. Pell and L. E. Pavne, 

Mathematika 7, 78-92 (1960). 
Reply to criticisms concerning Ireland contained in article, 

basic research in Europe, D. M. Gates, Science 138, No. 

3318, 1-9 (Aug. 1, 1958). 
Thermal voltage converters for accurate voltage measure- 
ments to 30 megacycles per second, F. L. Hermach and 

E. S. Wilhams, Commun. Electron. AIEE, No. 49, 200 

(July 1960). 
Modern theories of materials, C. Truesdeh, Trans. Soc. 

Rheology IV, 9 (1960). 
Use of the incoherent scatter technique to obtain ionospheric 

temperatures, T. E. VanZandt and K. L. Bowles, J. 

Geophys. Research 65, 2627 (Sept. 1960). 
Highlights of the consultative committee on international 

radio (CCIR) activities in the field of radio propagation, 

J. W. Herbstreit, 9th Plenarv Assembly of CCIR Proc. 

IRE 48, No. 1, 45 (Jan. 1960). 
Graphs for bivariate normal probabilities, M. Zelen and N. C. 

Severo, Ann. Math. Stat. 31, No. 3, 619 (Sept. 1960). 
Temperature stratification in a non-venting liquid hehum 

dewar, L. E. Scott, R. F. Robbins, and B. W. Birmingham, 

Proc. 1959 Cryogenic Engr. Conf. 5, (Sept. 2-4, 1959). 
A microwave impedance meter capable of high accuracy, 

R. W. Beatty, IRE Trans. Microwave Theory and Tech. 

MTT-8, No. 4 (July 1960). 
Variational treatment of electron-hydrogen atom elastic 

scattering, S. Geltman, Phys. Rev. 119, No. 4, 1283-1290 

(Aug. 1960). 
Relaxation processes in multistate systems, K. E. Shuler, 

Phys. Fluids 3, No. 4, 442 (July- Aug. 1959). 
Comment on models of the ionosphere above hmax F2, J. W. 

Wright, J. Geophys. Research 65, 2595 (Sept. 1960). 
Some magnetoionic phenomena of the Artie E-region, J. W. 

Wright, J. Atmospheric and Terrest. Phys. 18, 276 (Aug. 

1960). 
Tables of thermodynamic and transport properties of air, 
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gen, oxygen and steam, J. Hilsenrath, C. W. Beckett, 

W. S. Benedict, L. Fano, H. J. Hoge, J. F. Masi, R. L. 

Nuttall, Y. S. Touloukian, and H. W. Woolley, 478 pages 

(Pergamon Press, Oxford, London, New York, and Paris, 

1960). 
Methods of nuclear orientation, E. Ambler, Progress in 

Cryogenics 2, 235 (Heywood & Co., Ltd., London, England, 

1960). 
A conference on the propagation of ELF electromagnetic 

waves, J. R. Wait, Proc. IRE 48, 1648 (Sept. 1960). 
The strength of ten structural adhesives at temperatures 

down to minus 424 F., W. M. Frost, Proc. 1959 Cryogeni(? 

Engr. Conf. 5, 375 (Sept. 2-4, 1959). 
The Stokes flow about a spindle, W. H. Pell and L. E. Payne, 

Quart. Appl. Math. 18, 257 (1960). 
Possibihty of detecting ionospheric drifts from the occurrence 

of spread F echoes at low latitudes, R. W. Knecht, Nature 

Letter 187, 927 (Sept. 1960). 
A compilation and correlation of the PVT data of normal 

hydrogen from saturated liquid to 80 degree K., R. B. 

Steward and V. J. Johnson, Proc. 1959 Cryogenic Engr. 

Conf. 5, 548 (Sept. 2-4, 1959). 
Pressure dependence of rotationally perturbed lines in the 

ultraviolet band spectrum of CN., H. P. Broida and S. 

Golden, Can. J. Chem. 38, 1666 (1960). 
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molecules. IV. Solvent effect on internal rotarions, S. 

Lifson and I. Openheim, J. Chem. Phys. 33, No. 1, 109 

(July 1960). 
Improvements in radio propagation prediction service, W. B. 

Chadwick, Elec. Engr. 1-4 (Sept. 1960). 
Carbon resistance thermometrv with mixed dc and rf currents. 
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A. E. Boyd and E. E. Morris, Health Phys. 3, 321 (1960). 

Maser frequency stabihty, R. C. Mockler and J. A. Barnes, 
Proc. 13 Annual Frequency Control Symp., 583 (May 
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J. Am. Chem. Soc. 83, 5033 (1960). 
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dioxide powders; its dependency on porosity and grain size 
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The compound BaTiGeijO.), C. R. Robbins, J. Am. Ceram. 
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